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ABSTRACT 

We revisit the Blandford-Znajek process and solve the fundamental equation that governs the structure of the 
steady-state force-free magnetosphere around a Kerr black hole. The solution depends on the distributions of the 
magnetic field angular velocity co and the poloidal electric current I. These are not arbitrary. They are determined 
self-consistently by requiring that magnetic field lines cross smoothly the two singular surfaces of the problem: the 
inner “light surface” located inside the ergosphere and the outer “light surface” which is the generalization of the 
pulsar light cylinder. We find the solution for the simplest possible magnetic field configuration, the split monopole, 
through a numerical iterative relaxation method analogous to the one that yields the structure of the steady-state 
axisymmetric force-free pulsar magnetosphere. We obtain the rate of electromagnetic extraction of energy and 
confirm the results of Blandford and Znajek and of previous time-dependent simulations. Furthermore, we discuss 
the physical applicability of magnetic field configurations that do not cross both “light surfaces.” 
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1. INTRODUCTION 

The electromagnetic extraction of energy from a spinning 
black hole can be one of the most powerful and efficient engines 
in the universe. It was Blandford & Znajek (1977, hereafter 
BZ77) who first argued that a spinning black hole threaded 
by a sufficiently strong magnetic field and permeated by an 
electron-positron plasma generated from pair cascades estab- 
lishes a force-free magnetosphere, in direct analogy to the the- 
ory of axisymmetric pulsar magnetosphere developed a few 
years earlier by Goldreich & Julian (1969). They derived the 
fundamental equation that governs the structure of this mag- 
netosphere, the general relativistic force-free Grad-Shafranov 
equation, and obtained a second-order perturbative solution for 
small values of the Kerr parameter a and a scaling for the elec- 
tromagnetic power extracted from the rotating black hole. 

The structure, singular surfaces, and general properties of 
this equation for magnetic fields threading Schwarzschild and 
Kerr black holes were discussed recently by Uzdensky (2004, 
2005). He emphasized the need for a self-consistent determina- 
tion of the distributions of the magnetic field angular velocity 
and the poloidal electric current, something that early relaxation 
methods (Macdonald 1984) and later higher order perturbative 
solutions (Tanabe & Nagataki 2008) did not take into account. 
He also noted that this may be a very difficult task, and, indeed, 
to the best of our knowledge a solution of this equation in the 
general case of an open magnetic field configuration threading a 
Kerr black hole is still lacking. Nevertheless, the broader astro- 
physical interest in the Blandford-Znajek process has prompted 
in the past decade several researchers to seek and obtain the 
steady-state structure of the spinning black hole magneto- 
spheres through time-dependent general relativistic force-free 
numerical simulations (e.g., Komissarov & McKinney 2007; 
Palenzuela et al. 2011; Lyutikov & McKinney 2011). These 
works established the universal applicability of the electromag- 
netic extraction of the rotation energy of a Kerr black hole in 
powering active galactic nuclei, jets in X-ray binaries, and even 
gamma-ray bursts. 

The presence of this extra source of energy in addition 
to that of the accretion disk may account for some of the 


phenomenology. One should bear in mind, however, that the 
situation is more complicated. For instance, (1) almost 90% 
of extragalactic sources are radio quiet, i.e., they do not 
produce electromagnetic outflows in the form of radio jets, 
and yet they are believed to contain both constituents of 
the Blandford-Znajek process, namely a spinning black hole 
and a strong magnetic field threading it (e.g., Kukula 2003); 
(2) the formation and disruption of radio jets in black hole X-ray 
binaries over short timescales are at odds with this general notion 
(e.g., Belloni 2010); (3) there seems to exist strong observational 
evidence that the power in the radio jet is not directly related to 
black hole spin, contrary to the main result of BZ77 (Fender et al. 
2010; however, see also Narayan & McClintock 2012). Note, 
though, that this may be due to an erroneous determination of 
the black hole spin. 

The above justify an independent re-evaluation of the 
Blandford-Znajek process through the solution of the general 
relativistic force-free Grad-Shafranov equation. The special rel- 
ativistic form of that equation, the so-called pulsar equation, was 
only solved in 1999 by Contopoulos, Kazanas, and Fendt who 
showed that the magnetic field structure can only be obtained 
together with the determination of the unique electric current 
distribution that allows for a smooth and continuous structure 
to exist. In that respect, steady-state solutions offer a deeper 
insight into the physics of the problem than numerical time- 
dependent ones. It should be noted that, while BZ77 provided 
an estimate of the energy extraction efficiency of this process, 
they did not obtain an exact solution of the structure of the black 
hole magnetosphere. Also, at that time, the critical role of the 
singular surfaces and the distributions of the magnetospheric 
electric current and field line angular velocity had not been ap- 
preciated. Today, we revisit this problem with all the knowledge 
we carry from our 13 year long investigation of the force-free 
pulsar magnetosphere. 

2. THE GENERAL RELATIVISTIC PULSAR EQUATION 

In order to derive the fundamental equation that governs the 
steady-state structure of the force-free magnetosphere around 
a Kerr black hole, we follow closely the 3+1 formulation of 
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Thorne & Macdonald (1982) used by most researchers in the 
astrophysical community (e.g., Uzdensky 2005). We restrict 
our analysis to steady-state and axisymmetric spacetimes where 
(. . .)., = (. . ,)_0 = 0. In that case, the general four-dimensional 
spacetime geometry may be written in Boyer-Lindquist spher- 
ical coordinates x M = (f, r, 9, <p ) as 

ds 2 = g^dx^dx 11 

o , t A sin 2 9 „ , X , , 

= —a 2 dt 2 + {dtp - Q dtf + -dr 2 + Zd6 2 . (1) 

X A 

Here, a = (AX/ A) 1/2 and £2 = 2 aMr/A are the lapse function 
and angular velocity of “zero-angular momentum’' fiducial 
observers (ZAMOs), respectively, 

A = r 2 — 2Mr + a 2 , X = r 2 + a 2 cos 2 9, 

A = (r 2 + a 2 ) 2 — o 2 Asin 2 9, (2) 

where M is the black hole mass and a is its angular momentum 
(0 ^ a < M). Throughout this paper we adopt geometric units 
where G = c — 1 . Semicolon stands for covariant derivative and 
comma stands for partial derivative. Latin indices denote spatial 
components (1-3), Greek indices denote spacetime components 
(0-3), and denotes the spatial part of vectors. ZAMOs 
move with 4-velocity U u = ( I /a, 0. 0. il/a) orthogonal to 
hypersurfaces of constant t. The force-free magnetosphere of 
a spinning black hole is characterized by the electromagnetic 
energy-momentum tensor 



and the condition T.^ v = 0. Here, the rest mass and pres- 
sure contribution have been neglected. The electromagnetic 
field tensor F IUI is related to the electric and magnetic fields 
£''■ , B 1 ' measured by ZAMOs through F ilv — U tL E v — 
ljv E r + € ^pb x U p (e (jLvkp = [/x v Ap ] | det(g- M v ) | “ 1 /2 is the four- 
dimensional Levi-Civita tensor). Under these conditions, the 
fundamental equation that governs the steady-state structure of 
the force-free magnetosphere around a Kerr black hole becomes 

p e E + JxB = 0. (4) 

Here, p e and J are the electric charge and current densities, 
respectively. Equation (4) is supplemented by Maxwell’s equa- 
tions of electrodynamics 


V • B = 0 

V ■ E = 47 Xp e 
V x {a B) = Ana] 


and the electric and magnetic vector fields can be expressed in 
terms of three scalar functions, v P(r, 8), coQV), and 7(T / ) as 


B(r,6)= — (9) 

VAsin0 [ « 

E(r,6)= (10) 

aVx 

where to is the angular velocity of the magnetic field lines, 
I is the poloidal electric current flowing through the circular 
loop r — const., 9 = const., and is equal to (2;r) _1 times 
the total magnetic flux enclosed in that loop. Note that the 
electric field changes sign close to the horizon with respect 
to its sign at large distances. As explained in BZ77, a rotating 
observer (ZAMO) will in general see a Poynting flux of energy 
entering the horizon, but he will also see a sufficiently strong 
flux of angular momentum leaving the horizon. That ensures that 
energy is extracted from the black hole. The poloidal component 
of Equation (4) then yields the general relativistic force-free 
Grad-Shafranov equation: 


1 

1 
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to 2 A sin 2 9 AMator sin 2 9 
1 + - 
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cos 9 
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to' sin 2 9 


-(coA - 2a Mr) I v Py + - V U 


= II' 

A 

( 11 ) 


(Equation (3.14) of BZ77 rewritten in our notation). Henceforth, 
primes will denote differentiation with respect to V P. One sees 
directly that if we set a = 0 and M — 0 in Equation (11), we 
obtain 


TL 


+ — 'P, 00 + 


2M J , 


1 COS0 


ML 


r- 

2M J , 


x [1 — trr 2 sin 2 9 ] 


2to 2 cos 9 sin d'A’ y — toto'r 2 sin 2 9 


x j H' 2 + -T 2 ) = -All', 


( 12 ) 


V x (uE) = 0. 


Here, 


(5) which is the well-known pulsar equation (Scharlemann & 
Wagoner 1973). The zeroing of the term multiplying the second- 
order derivatives in Equation (11), 


V • A - A !■ . (V x Ay = [ijk\\d S t(g lm )\- l ' 2 A t j , (6) 

A - B = g'-i AjBj, (A x By = [ijk]\det(g, m )\- l ' 2 AjB k . (7) 

For several applications in astrophysics, perfect (infinite) con- 
ductivity is a valid approximation. In this case, 

E ■ B = 0, (8) 


to 2 A sin 2 9 AMator sin 2 9 2 Mr 

1 + — — = 0 , (13) 


yields the singular surfaces of the problem. We will henceforth 
call the singular surfaces “light surfaces” (LS). 4 When we 

4 These are none other than the “velocity-of-light surfaces” of Macdonald 
(1984) where the speed of a particle moving purely toroidally with angular 
velocity co equals c. 


2 


The Astrophysical Journal, 765:113 (6pp), 2013 March 10 


Contopoulos, Kazanas, & Papadopoulos 


set M = a = 0, it yields the standard pulsar light cylinder 
rsin0 = c/co. In general, the shape of the outer LS is only 
asymptotically cylindrical as 9 -* 0 (see Figure 2), and an inner 
LS appears inside the ergosphere. We will henceforth use the 
notation “light cylinder” (LC) and “outer LS” interchangeably. 
It is interesting to note that the outer boundary of the ergosphere 
corresponds to the solution of the singularity condition for 
co = 0, whereas the inner boundary (the event horizon at 
r = /- BH = M + V M 2 — a 2 ) corresponds to the solution of 
the singularity condition for co = Qbh = a/(r BH + a2 )> where 
£2 bh is the angular velocity of the black hole. It is also interesting 
to note that the natural “radiation condition” at infinity (energy 
must flow outward along all field lines) requires that 

£2 B h (14) 

(BZ77), and therefore indeed the inner LS lies inside the 
ergosphere. 

Both Equations (11) and (12) contain the two functions, 
o>( y Y) and /(M'), which must be determined by the physics of 
the problem. In the case of an axisymmetric spinning neutron 
star, co is usually taken to be equal to the neutron-star angular 
velocity In pulsars, / (*F) is self-consistently determined 

through an iterative numerical technique that implements a 
smooth crossing of the relativistic Alfven surface, which is 
the usual light cylinder where rsind = c/co (Contopoulos 
et al. 1999; Timokhin 2006). In the case of a spinning black 
hole, the situation is qualitatively similar but more complicated. 
Contrary to a neutron star, the black hole does not have a 
solid surface, and therefore it cannot impose any restriction 
on co other than the “radiation condition” Equation (14). The 
only natural restriction that is imposed by the physical problem 
is that magnetic field lines must be smooth and continuous 
everywhere. co( y Y) must therefore be determined together with 
IC ¥ ) through the condition of smooth crossing of both LS of 
the problem: the inner one inside the ergosphere and the outer 
one at large distances. In the next section, we will see that 
this can be achieved through an iterative numerical technique 
analogous to the one employed in the pulsar magnetosphere. As 
remarked previously, iterating with respect to two functions 
simultaneously is indeed a very difficult task. This is why 
Uzdensky (2005) opted to consider only field lines that connect 
the black hole horizon to the surrounding accretion disk where 
<w( v F) is determined by the Keplerian angular velocity of the 
disk. 

There exists one more interesting complication. Equation (14) 
applies only to open field lines that cross the event horizon. In 
the case of a star, field lines that cross the stellar surface rotate 
with the angular velocity of the star, but they do not cross an 
event horizon (either because no event horizon forms or because 
they simply avoid the horizon), and the above restriction does 
not apply. In the present case, we are interested in studying the 
fundamentals of the Blandford-Znajek process, and therefore, 
we will consider a black hole “as clean as possible.” Obviously, 
a magnetized astrophysical black hole is not isolated. In that 
respect, we would like here to study an astrophysical black 
hole with the minimum number of extra elements, namely 
a surrounding thin disk of matter to hold the magnetic field 
through the necessary electric currents and charges. Charges 
and currents will be produced through pair production in the 


5 Note that in the presence of particle acceleration magnetospheric “gaps,” 
this is not 100% exact (Ruderman & Sutherland 1975; Contopoulos 2005). In 
particular, in old pulsars near their death line cj <$C Qns • 


black hole magnetosphere. We will ignore magnetic field lines 
that do not cross both the inner and outer LS because we will 
have no physical way to determine their co and the poloidal 
electric current I contained inside them. We will also ignore 
magnetic field lines that cross the equator outside the black hole 
event horizon because such field lines need a source of poloidal 
electric current on the equator (if they rotate, they will need to 
cross at least one LS where the condition of smooth crossing 
will in general require a non-zero electric current to flow along 
that line). 

Therefore, in the present work we will only study black hole 
magnetospheres where all magnetic field lines reach the event 
horizon, and all magnetic field lines are open. As argued in 
Lyutikov (2012), such configuration is the most natural end 
stage of stellar collapse. Note once again that it is sustained 
by electric currents in an equatorial current sheet. The effect of 
a surrounding disk of a certain height possibly threaded by the 
return magnetic field (as in the scenario predicted by the Cosmic 
Battery; Contopoulos & Kazanas 1998) will be considered in a 
future work. 


3. THE NUMERICAL ITERATIVE SCHEME 


We will here describe the technical numerical details that led 
to the solution of Equation (11). First, we change the radial 
variable from r to R(r) = r/(r + M). Our numerical integration 
extends from R mm = R(r B h) (the event horizon) to R max = 1 
(radial infinity). The 6 coordinate extends from 0 tmn = 0 (the 
axis of symmetry) to 0 max = jt/ 2 (the equatorial plane). We 
implemented a 256 x 64 numerical grid uniform in R and 9. 
Note that this grid has a very high resolution in r around the 
black hole event horizon where the inner LS lies, but not as high 
around the outer LS, and in particular at low 9 and a values. 

Second, we specify boundary conditions on the axis of 
symmetry, the horizon, the equatorial plane, and infinity. V P(0 = 
0) = 0. TO'bh, 9) is determined by the integration over 9 of the 
so-called Znajek regularization condition 


im = 


Mr Bli sin 9 

r BH + fl2 cos2 $ 


(^BH 


- coy¥, e 


(15) 


(Znajek 1977). Note that this condition is updated as the 
distribution /( V P) is updated along with the numerical solution 
for v P(r, 9). The equatorial boundary condition is interesting. 
As we discussed in the previous section, the equatorial region 
may contain an accretion disk, with or without its own magnetic 
field. Before investigating such a complex astrophysical system, 
we would like to understand first the simplest case, that of a 
magnetized black hole where the electric currents supporting 
its magnetic field are distributed on a thin equatorial current 
sheet. In that case, the equatorial boundary condition becomes 
v P(r, 9 = jr/2) = TVix- The outer boundary condition is also 
not important since the outer boundary has practically been 
moved to infinity. In practice, we chose '¥. I (R = 1) = 0, but any 
other boundary condition yields the same results. We initialize 
our numerical grid with a split monopole field configuration 


l F(r, 9) = ‘FsmCl 9) = TW1 - cos 9) , (16) 


with 


aifP) = 0.5£2 bh 

and 

I OP) = Ism m = -0.25 Q bh ‘F[ 2 - (‘P/‘P max )]. (17) 
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The reader can directly check that Equations ( 1 6) and (17) are an 
exact analytical solution of the pulsar equation (Equation (12); 
Michel 1982). Note the negative sign in Equation (17). For the 
particular magnetic field configuration where B' on the axis is 
aligned with the black hole spin direction, this corresponds to 
an outflow of electrons. Different initial configurations yield the 
same final solution. We discretize all physical quantities on our 
256 x 64 ( R , 0) grid and we update V F( R , 0) through simulta- 
neous overrelaxation with Chebyshev acceleration (subroutine 
SOR from Numerical Recipes; Press et al. 1986). 

Third, we update the distributions of o>( l F) and IC¥) as 
follows: at each latitude 9, we check where the singularity 
condition (Equation (13)) is satisfied in r. At each such radial 
position, we extrapolate l F inward from larger r (M'(r + , 9)) and 
outward from smaller r ( v F(r~, ())). In general, v P(r + , 0) and 
v F(r~, 9) differ. Then, at the inner LS we implement 

^newCPnew) = ^old^new) - 0.1[T'(r + , 9) - 0)] , (18) 

whereas at the outer LS we implement 

/newOPnew) = /oldC^new) + 0.1[TV, 9) - T(r", 0)] 

WnewOPnew) = ^oldCPnew) - 0.1[‘P(r + , 0) - TV", #)], (19) 

where 

TV e w = O.5[‘F(r + ,0) + T'(r-,0)] (20) 

at each LS. The reasoning here is that we impose weighted 
corrections on <y( l F) and / PF) based on the non-smoothness of 
the *P(r, 0) distribution along all grid points inside and outside 
the two LS, and not through the regularization conditions as we 
did in Contopoulos et al. ( 1 999). This is a very general procedure 
that may be applied to any similar singular equation. 

Our numerical scheme does not converge easily, and its 
parameters need to be empirically adjusted by following the 
progress of the iteration. The signs of the weighting coefficients 
are determined by trial and error: when we choose the wrong 
signs, the mismatches | x F(r + , 0) — 0) | increase and the 

*P(r, 0) distribution is very quickly destroyed. Their magnitudes 
are also obtained empirically: too large values (around unity) 
lead to overcorrections and instability and too small values 
(below 0.01) do not correct fast enough for the iteration to 
converge. Furthermore, in order to facilitate convergence, at 
the inner LS we update only a>, whereas at the outer LS we 
update both I and co at every other relaxation iteration for X P. 
Also, in order to avoid numerical instabilities, we smooth out 
the distributions of l and co every 50 relaxation iterations for Mb 
We found that more frequent smoothing inhibits convergence. 
Our goal is to minimize the residuals in the discretized form of 
Equation (11) together with the non-smoothness in MTr, 0), but 
our method fails for poor grid resolutions at the inner and/or 
the outer LS (see Discussion section). 

The results of our numerical integration are shown in Figure 1 . 
We show here the distributions of &>( l F) (top panel), 2|/( V P)| 
(center panel), and Poynting flux at large distances 2|/(T / )|a>(T / ) 
(bottom panel) after 6000 relaxation steps for X ¥{R,9), for 
six values of the black hole spin parameter a. As is the 
case in pulsars, the magneto spheric electric current is non- 
zero at T = v F lnax , and the global electric circuit closes 
through an equatorial current sheet. Note that we have no way 
to update co along the axis, and therefore we have chosen 
(oBY = 0) = 0.5£2bh- 6 On the equator we implemented 

6 We have also ran our simulations with co'QV = 0) = 0 and the solutions 
were similar, only the convergence of a> became more unstable. 



0 0.2 0.4 0.6 0.8 1 




Figure 1 . Distributions of &i(4 / ) normalized to Qbh (top panel), 2|/('F)| 
normalized to flnn'Lmix (center panel), and Poynting flux at large distances 
2|/('P)|a)('t / ) normalized to 12/ r[ 'l\ Tlax (bottom panel) for various values of the 
black hole spin parameter: a/M = 0.7 (solid), 0.8 (short dashed), 0.9 (dotted), 
0.99 (small dots), 0.999 (dash dotted), and 0.9999 (short dash dotted). All 
three distributions are close to the analytical split monopole expressions of 
Equation (17) (thick dots). 


4 


The Astrophysical Journal, 765:113 (6pp), 2013 March 10 


Contopoulos, Kazanas, & Papadopoulos 




Figure 2. Poloidal magnetic field lines (lines of constant 'P in 20 equal subdivisions of 'P max ) for a/M = 0.9999 close to the event horizon (left panel) and farther out 
(right panel). The thick lines are the event horizon, the inner LS, the outer boundary of the ergosphere, and the outer LS (from the center outward, respectively). 


co'Q V = M' nlax ) = /'(M' = M' rna J = 0. The reader should 
check the qualitative agreement with the results of the time- 
dependent force-free general relativistic simulations performed 
by Tchekhovskoy et al. (2010, their Figure 4). 

In Figure 2 we show the characteristic magnetic field con- 
figuration for a/M — 0.9999 near the black hole horizon (left 
panel) and at larger distances (right panel). We also plot the 
inner LS inside the ergosphere, and the outer LS which be- 
comes asymptotically cylindrical as 6 — > 0. Note that magnetic 
field lines are monopolar at large distances (beyond a few times 
the radius of the event horizon), and they bend upward toward 
the axis as they approach the horizon. This effect is also seen 
in the numerical simulations of Lyutikov & McKinney (2011) 
for a/M = 0.99. It is not discemable in the field configura- 
tions of Figure 1 of Tchekhovskoy et al. (2010) because the 
scale is not adequate, but is implied in their Figure 7. Note that 
ftj(T') is equal to 0.5£2bh to within 10%, and that 7(*P) is very 
close to the split monopole expression /sm( x P) also to within 
10% (Equations (16) and (17)). Therefore, the Znajek condition 
(Equation (15)) can be rewritten using Equation (9) as 


B 1 


VI/ 

^ 1 max 

lr=r BH ~ 2 Mr BH 

A max 


(2 M) 2 


Va 

i + 


(a/M) 1 


[1 + i/l — {a/M) 2 ] 1 


cos 2 0 


( 21 ) 


Finally, we calculate the total rate of electromagnetic energy 
extraction (from both hemispheres) as measured at infinity 
through Equation (4.1 1) of BZ77 as 


/■'{'max 

8= -2x2 tt / 2 1('¥)coQ¥)d'¥ 

J'V=0 

'A 2 it/2 


— kSsM — k ^ ^BH^max- 


( 22 ) 


We have here normalized our results to the split monopole 
value with cusm = 0.5£2bh- k & 0.8 for a/M — 0.9999 and 
is practically indistinguishable from unity for a/M 0.9 (see 
the bottom panel of Figure 1), thus confirming the main result of 
BZ77, namely that, in an open magnetic field configuration with 
an infinitely thin surrounding equatorial disk, electromagnetic 
energy is extracted at a rate proportional to both Q| h and T ,2 lax . 


4. DISCUSSION 

We believe that the study of the steady-state open force- 
free black hole magnetosphere offers a deeper insight into 
the physics of the problem than time-dependent numerical 
simulations. In analogy to the axisymmetric force-free pulsar 
magnetosphere, we argued that steady-state solutions can only 
be obtained through a self-consistent determination of the 
distributions of both the magnetic field angular velocity and the 
magnetospheric electric current cu(T / ) and 7( X T / ), respectively. 
These can only he determined when open field lines cross both 
singular LS: the inner LS inside the ergosphere outside the black 
hole horizon and the outer LS which is a deformed light cylinder. 
We developed an iterative numerical procedure analogous to the 
one we implemented in the study of the pulsar magnetosphere. 
Iterating with respect to two functions simultaneously is a very 
difficult task and our method may certainly be improved in the 
future. However, we are confident enough to suggest that the 
distributions are unique for the particular boundary conditions 
of the problem (all magnetic field lines thread the black hole 
horizon, all are open, and the surrounding disk is infinitely thin) 
since different initial conditions for ¥, a>, and 7 evolve toward 
the same final solution (for each value of a that is). Our solutions 
may serve as test cases for time-dependent numerical codes. We 
plan to apply our method to the study of various astrophysically 
interesting situations where the surrounding disk is not infinitely 
thin but expands vertically with distance. 

All field lines that cross the black hole horizon are bound to 
cross the inner LS. This is not the case, though, for the outer 
LS. In fact, several numerical simulations exist in the literature 
where the black hole is threaded by a vertical uniform magnetic 
field, and yet they converge to certain distributions of o>( V P) and 
7(T / ) (e.g., Komissarov & McKinney 2007; Palenzuela et al. 
2011). We expect that in that case, the solution is mathematically 
indeterminate, and the solutions shown depend on the particular 
choice of outer numerical boundary conditions. For example, 
the lateral boundary conditions in Palenzuela et al. (2011) 
are not perfectly absorbing but reflect outgoing waves as is 
manifested in the oscillations seen in their Figure 2. We plan to 
study such magnetic field configurations in a future work and 
check whether we can indeed freely specify one of the two free 
functions, &>( V F) or / ( V F). We will thus check whether the results 
of previous time-dependent numerical simulations are unique 
or not. 
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We would like to caution the reader about certain numerical 
complications that arise in all numerical simulations of the 
spinning black hole magnetosphere, time dependent and steady 
state. Numerical simulations that do not adequately resolve 
both the inner and outer LS should not be trusted. We faced 
a similar problem in our study of the pulsar magnetosphere. 
In that case, there exists no inner LS, and in order to study 
the transition through the outer LS (the light cylinder) we 
were forced to implement our inner boundary condition at 
a radial distance more than a hundred times larger than the 
actual neutron-star radius (Spitkovsky 2006; McKinney 2006; 
Contopoulos & Kalapotharakos 2010). This, however, posed 
no problem since the angular velocity of the field lines is set 
by the rotation of the pulsar. The black hole problem is more 
complicated because the black hole surface is not infinitely 
conducting, and the numerical code must adequately resolve 
both the inner boundary and the light cylinder. Note that our 
numerical method works best as a / M approaches unity since in 
that case the inner LS is well separated from the event horizon, 
and the outer LS is at its closest approach to the origin of 
the numerical grid where its radial resolution is the highest. 
In that case, a numerical resolution of 256 radial grid points 
corresponds to 34 grid cells between the horizon and the inner 
LS, and 94 grid cells between the two LS along the equator. 
When a/M 2; 0.9, our results are rather resolution independent. 
On the other hand, as a/M approaches zero, the ergosphere 
shrinks in width and the inner LS approaches so close to the 
event horizon that our numerical grid does not have enough 
radial resolution for the iterative procedure that yields cuRE) 
to properly work. 7 At the same time, as a/M decreases so 
does £2 bh, pushing the outer LS so far out that our radially 
expanding numerical grid has very low radial resolution at the 
position of the outer LS for the iterative procedure to properly 
converge. Similarly, the resolution of our radially expanding 
grid on the light cylinder becomes worse as 6 — > 0, and this 
too complicates the convergence of our iterative method. Other 
numerical approaches certainly face similar complications at 
short and/or at large distances. Our method may be improved 
by a reformulation in radially expanding cylindrical coordinates, 
in Kerr-Schild coordinates, or in standard spherical coordinates 
with a numerical integration outer radial boundary. 

In summary, the force-free spinning black hole magneto- 
sphere is not very different from the force-free axisymmetric 
pulsar magnetosphere. In fact, beyond a distance on the order of 
a few times the black hole radius, the two problems are almost 
identical. An important corollary of this similarity is that a spin- 
ning black hole surrounded by a thin equatorial plasma disk 
does not generate a collimated outflow, only an uncollimated 
wind (as in Contopoulos et al. 1999). Thus, the problem of the 
black hole jet formation remains open. A possible resolution 
must involve a surrounding thick disk, or a surrounding mag- 
netic disk wind that gradually collimate the initially monopolar 
black hole wind into a collimated jet outflow, as seen in the 


numerical simulations of Tchekhovskoy et al. (2010). It is in- 
teresting to note here, again similarly to pulsars, that beyond 
the light cylinder the Lorentz factor T in the monopolar wind 
increases linearly with cylindrical distance as 

r sin 0 rO BH sin 6 

T oc ^ — — (23) 

c/to 2c 

(Contopoulos & Kazanas 2002), possibly up to the fast 
magnetosonic point of the outflow. Therefore, according to 
Equation (23), collimated field/flow lines that do not extend 
as far from the axis as uncollimated ones are expected to reach 
lower Lorentz factors. Of course, the final Lorentz factor of such 
outflows will also depend on the mass loading of a given mag- 
netic field lines, but such considerations are beyond the scope 
of the present work. We plan to continue our investigation of 
the interrelation between the black hole and disk outflows in the 
future. 
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